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Abstract 

In this paper, the author present a reliable symbolic computational algorithm for inverting a general 
comrade matrix by using parallel computing along with recursion. The computational cost of our 
algorithm is 0(n 2 ). The algorithm is implementable to the Computer Algebra System (CAS) such as 
MAPLE, MATLAB and MATHEMATICA. Three examples are presented for the sake of illustration. 
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1 Introduction 

The n x n general comrade matrix, denoted by C, takes the form 
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The comrade matrix is a generalization of the companion matrix and is associated with a polynomial ex- 
pressed as a linear combination of an arbitrary orthogonal basis. This matrix appears frequently in many 
areas of science and engineering, for example in linear multivariable systems theory [1], Computing the 
Greatest Common Divisor of Polynomials [2] and division of generalized polynomials [3]. The solution of 
comrade linear system has been investigated by many authors (see for instance, [4-6]). Finding the inverse 
of comrade matrix is usually required to solve this linear system. The motivation of the current paper is 
to establish efficient algorithm for inverting companion matrix of the form (1.1). 

The paper is organized as follows. In Section 2, new symbolic computational algorithm, that will not 
break, is constructed. In Section 3, three illustrative examples are given. Conclusions of the work are 
presented in Section 4. 
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2 Main Results 



In this section we shall focus on the construction of new symbolic computational algorithm for computing 
the determinant and the inverse of a general comrade matrix. Firstly, we begin with computing the LU 
factorization of the matrix C. It is as in the following: 
Let 

C = LU 
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The elements in the matrices L and U in (2.2) and (2.3) satisfy: 
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We also have: 



Xi = { —(a n -i + i — Ui^\Xi-\) 
^77(7™ - a n - 2 x n - 2 ) 

Det(C) 



if i 

if i = 2,3,..., n- 1 
if i = n, 

if i = l 

if i = 2, 3, n - 
if i = n — 1. 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



At this point it is convenient to formulate our first result. It is a symbolic algorithm for computing the 
determinant of a comrade matrix C of the form (1.1). 

Algorithm 2.1 To compute Det(C) for the comrade matrix C in (1.1), we may proceed as follows: 
INPUT order of the matrix n and the components on, i = 1,2, ...,n — 1, Pi, i = 1,2, ...,n, 

7i, i — 2, 3, n, and ai, i = 3,4, n. 
OUTPUT The determinant of comrade matrix C. 

Step 1: Set = P\. If fix — then /ii = t(t is just a symbolic name) end if. Set x\ 
Step 2: Fori = 2,3,...,n-2 

Compute \ii — (3i — -p—^oii-i, if fii = then [ii = t, 

Compute x t = ^-(a n _ i+ i - cti-iXi-i), 
Step 3: Set ^„_i = /3„_i - ?f^-a n - 2 , «/m«-i = then yu n _i 



mi 



t end if, 



Set x n -i = j^(jn - x n - 2 a n -2) ■ 



Set fj, n = Pn ~ Otn-lXn-\, if Hn = then fl r , 

Step 4-' Compute Det(C') — (^Yl7=i 



t end if. 



The symbolic Algorithm 2.1 will be referred to as DETSGCM algorithm. The computational cost of 
DETSGCM algorithm is 7n - 10 operations. The new algorithm DETSGCM is very useful to check 
the nonsingularity of the matrix C 



Now, when the matrix C is nonsingular, its inversion is computed as follows: 

Let 

C' 1 = [Sij]i<i^< n - [Col u Col 2 ,...,Col n }. (2.7) 
where Col m denotes m t h column of C _1 , m — 1,2, ...,n. 

Since the Doolittle LU factorization of the matrix C in (1.1) is always possible then we can use parallel 
computations to get the elements of the last two columns Cok — (Si^, 62,,, S, h i) T , i = n and n — 1 of 
C- 1 as follows [7]: 

Solving in parallel the standard linear systems whose coefficient matrix L is given by (2.2) 
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we get 
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Hence, solving the following standard linear systems whose coefficient matrix U is given by (2.3) 
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gives the two columns Col^ i = n and n — 1 in the forms: 
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aiSi +^~\ i = „-2,n-3,...,l. (2.15) 

Mi 

Using equations (2.11)-(2.15) with the fact C~ 1 C = /„ where /„ is the n x n identity matrix, the 
elements in the remaining (n - 2) columns of C _1 may be obtained recursively using are obtained by using 

Col n -2 = —^—{En-x - /?„_iCo/„_i - 7„CoZ„), (2.16) 

a-n-2 

Colj = — (E j+1 - f3 j+1 Col J+1 - j j+2 Col j+2 - a n -jCol n ), j = n - 3,n- 4, 1. (2.17) 

Ctj 

Here E r = (Si r , S 2r , ■ S nr ) T , r = 1,2, n, where #i r is the Kronecker symbol. 

Now we formulate a second result. It is a symbolic computational algorithm to compute the inverse of 
a general comrade matrix of the form (1.1) when it exists. 



Algorithm 2.2 To find the n x n inverse matrix of the general comrade matrix C in (1.1) by using the 
relations (2.11)-(2.17). 

INPUT Order of the matrix n and the components a i} i = 1,2, n — 1, ft, i = 1, 2, n, 

7i, i = 2, 3, n, and (Zj, i = 3, 4, n. 
OUTPUT Inverse comrade matrix C~ x . 

Step 1: Use the DETSGCM algorithm to check the nonsingularity of the matrix C. If the matrix C is 
singular 

then OUTPUT (The matrix C is singular), Stop. 
Step 2: If ^ = for any i = 1,2, n, set pn = t (t is just a symbolic name). 
Step 3: If a.i = for any i — 1,2, n — 2, set a-i — t. 

Step 4-' For i = 1,2, n, compute and simplify the components Si >n and S'^n-i of the columns Colj, 

j = n, and n — 1, respectively, by using (2.11)-(2.15). 
Step 5: For i = 1,2, ...,n, compute and simplify the components Si >n -2 by using (2.16). 
Step 6: For j = n — 3, n — 4, 1, do 
For i = 1,2, n, do 

Compute and simplify the components Sij by using (2.17). 
End do 
End do 

Step 7: Substitute the actual value t = in all expressions to obtain the elements, Sij, i,j = 1, 2, ...,n. 

The symbolic Algorithm 2.2 will be referred to as SGCMINV. The computational cost of SGCMINV 
algorithm is 7n 2 — 5n — 11 operations. In [8], the author presented recurrence relations for the rows of an 
inverse comrade matrix but he supposed that on ^ for i = 1, 2, n — 1 and the first row of an inverse 
comrade matrix is known. The computational cost of this method is 0(n 3 ) operations. On the other hand, 
if we set eij = 0, i = 3, 4, n, the Algorithm 2.3 in [7] will be special case of the SGCMINV algorithm. 



3 ILLUSTRATIVE EXAMPLES 

In this section we give three examples for the sake of illustration. 
Example 3.1. Consider the 7x7 matrices C given by 



C = 



5 





1 

\ 
I 5 

3 



_ 1 

3 





3 

\ 

3 







1 

3 

3 



o \ 




1 

3 

3 



A 



by applying the SGCMINV algorithm. 
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(Steps 2-7). 



Example 3.2. Consider the 4x4 matrices C given by 
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by applying the SGCMINV algorithm. 
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(Steps 2-7). 



Example 3.3. We consider the following n x n comrade matrix in order to demonstrate the efficiency of 
SGCMINV algorithm. 
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We used SGCMINV algorithm to compute the inverse of comrade matrix C. Results are given in the 
next table in which e = \ \C~^ act 
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50 


100 


500 


: '- 1 Cexact ^SGCMINV 1 loo 


1.1631 x 10" 9 


1.1215 x 10~ 9 


1.6078 x 10~ 9 


CPU time(s) 


0.421 (using inverse function in Maple 13.0) 


3.448 


336.338 


0.109(Using our Algorithm) 


0.609 


33.899 



4 CONCLUSIONS 



In this work new symbolic computational algorithms have been developed for computing the determinant 
and inverse of a general comrade matrix. The algorithms are reliable, computationally efficient and remove 
the cases where the numeric algorithms are fail. 
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